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Abstract 

In this paper we explore how far the post-Newtonian theory, goes in overcoming the 
difficulties associated with anisotropic homogeneous cosmologies in the Newtonian approx- 
imation. It will be shown that, unlike in the Newtonian case, the cosmological equations 
of the post-Newtonian approximation are much more in the spirit of general relativity with 
regard to the nine Bianchi types and issues of singularities. 

The situations of vanishing rotation and vanishing shear are treated separately. The 
homogeneous Bianchi I model is considered as an example of a rotation-free cosmology with 
anisotropy. It is found in the Newtonian approximation that there are arbitrary functions 
that need to be given for all time if the initial value problem is to be well-posed, while 
in the post-Newtonian case there is no such need. For the general case of a perfect fluid 
only the post-Newtonian theory can satisfactorily describe the effects of pressure. This is 
in accordance with findings in [Q where the post-Newtonian approximation was applied to 
homogeneous cosmologies. 

For a shear-free anisotropic homogeneous cosmology the Newtonian theory of Heckmann 
and Schiicking, |^ is explored. Comparisons with its relativistic and post-Newtonian coun- 
terparts are made. In the Newtonian theory solutions exist to which there are no analogues 
in general relativity. The post-Newtonian approximation may provide a way out. 



1 



1 Introduction 



Due to its simplicity, the Newtonian approximation of cosmology is preferable to general rela- 
tivity and is typically used wherever possible. However, in a recent paper Q, it is shown that 
the Newtonian approximation has some difficulties. To begin with the Newtonian approxima- 
tion is incomplete in the sense that the Bianchi identities cannot be obtained from the field 
equations. It is generally accepted that Newtonian theory is a specialization of the linearized 
approximation of general relativity cut off at the c~^-level. However, information is lost in this 
cut-off procedure, and it is shown that cutting off at the c~^-level results in a theory in which 
the Bianchi identities are obtainable from the field equations, Moreover, the Newtonian 
theory has no initial value problem when applied to cosmology because there are no bound- 
ary conditions. Therefore, the Poisson equation has no unique solution. The post-Newtonian 
theory, on the other hand, is a closed, consistent and well-posed theory. The post-Newtonian 
approximation used here is obtained by expanding up to order and reformulating the lin- 
earized field equations as wavelike equations. A specialization of the resulting equations leads 
to a formally well-posed c~^-theory. This c~^-approximation is the first order approximation to 
give us consistency and formal well-posedness and, as we shall see, suffices as a "Newtonian" 
approximation to cosmology. Going to higher order simply gives increasingly more accurate 
approximations of general relativity. 

In 0] it was shown that the post-Newtonian approximation, in the context of the k = 
Friedmann-Robertson- Walker cosmology (FRW), is able to fully reproduce the results of its 
general relativistic counterpart, whilst the Newtonian theory is not. The Newtonian approxi- 
mation of cosmology does not allow the pressure to enter into the dynamics. Hence, changing 
the equation of state, does not result in different solutions for the density and expansion. In 
fact, the only solution that agrees with general relativity is that for the matter dominated uni- 
verse. On the other hand, in the post-Newtonian approximation pressure becomes dynamic 
through an extra function of time, A{t). Hence, the post-Newtonian theory is able to produce 
differing solutions for alternative equations of state in accordance with general relativity. In 
the case where A{t) goes to zero the pressure is forced to vanish and the higher order terms of 
the post-Newtonian approximation vanish as well, giving exactly the same solution for dust as 
the Newtonian theory does. Thus, Newtonian theory for homogeneous and isotropic cosmology 
should only be used for the special case of dust, where the two theories, the Newtonian and 
post-Newtonian, coincide. 

In this paper we explore how far the post-Newtonian approximation goes in overcoming 
the difficulties of the Newtonian theory when applied to the more general case of anisotropic 
homogeneous cosmology. There is a great deal of theoretical and observational evidence to 
support anisotropy in the universe, Any theory which attempts to approximate general 
relativity should yield similar results to the fully relativistic theory. The FRW models are 
isotropic and homogeneous. Although they fit most of the observed features of the universe 
they do not fit all. Different geometries may be worth considering at earlier times, for example, 
near the initial singularity. There might exist a general rotation of matter in the universe of 
comparable magnitude to the general expansion that we can not detect. Shear might provide 
a deviation away from isotropy since extra-galactic objects might be observed by the galactic 
equator Q. All homogeneous cosmologies fall into two classes: the Bianchi models, which are 
those for which the isometry group admits a 3-dimensional simply transitive subgroup, and 
those for which the isometry group is neither simply transitive, nor admits a simply transitive 
subgroup - the Kantowski-Sachs models. There exist nine Bianchi types and, correspondingly, 
nine Bianchi cosmologies, each class has subclasses with extra symmetries. The Bianchi types 
are in general anisotropic, so they do not have all spatial directions at a point being equivalent. 
Instead, there are now accelerations, and anisotropic observers will no longer be orthogonal to the 
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surfaces of homogeneity. Not all Bianchi types are suitable for describing the universe. They were 
first considered as cosmological models by Taub, |11| then shortly followed by Raychaudhuri, 
Q and Heckmann and Schiicking, [Q]. These models contain the FRW models as a special case. 
Thus, they are considered the simplest and the most likely cosmological candidates, in particular 
Bianchi I, Bianchi IX and the LRS types 7o and 7^. 

Bianchi I has a different cosmic scale factor for each direction and only becomes isotropic in 
the case where the scale factors all become equal. In this case the Bianchi I model becomes the 
A; = FRW cosmology. For nonzero shear these rotation-free anisotropic models are called the 
Heckmann- Schiicking solutions, [^.Q In the vacuum case they become the Kasner solution. 

Unfortunately, there are solutions in the case of anisotropy to which there are no general 
relativistic analogues. Theorems of Ellis, tell us that shear-free perfect fluids have either 
vanishing expansion or vanishing rotation. The case of vanishing rotation just yields the FRW 
model which is explored in Q. The case of vanishing expansion, however, is more problematic 
since in Newtonian cosmology it allows for singularity-free solutions. This in in contradiction 
with theorems of Hawking which state that solutions of general relativity must have singularities. 
The post-Newtonian approximation provides a way out of this difficulty. 

We begin in section 2.1 with a review of homogeneous Newtonian cosmology. Then in 2.2 
we consider post-Newtonian approximation to homogeneous anisotropic cosmology. We find 
that the Newtonian theory is in general under-determined and unable to fully reproduce the 
nine Bianchi types of general relativity. But the post-Newtonian equations produce a set of 
cosmological equations which are more in the spirit of the Bianchi identities. 

In section 3 we consider the case where the rotation vanishes separately. Bianchi I is taken 
as an example of a rotation-free anisotropic and homogeneous cosmology, and explored in the 
Newtonian and post-Newtonian contexts. In section 4 we consider the shear-free case. We end 
in section 5 with a conclusion. 



^Not to be confused with the Newtonian Heckmann-Schiicking cosmologies. 
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2 Newtonian Theory and the Post-Newtonian Approximation 
in Homogeneous Anisotropic Cosmologies 

2.1 Newtonian Homogeneous Anisotropic Cosmology 

In the following we will consider the Newtonian cosmological model of Heckmann and Schiicking, 
m. For a homogeneous cosmology the density, p, and the pressure, p, are functions of time only. 
The velocity field, Vi, is the same relative to all observers and such that Vi = Vij{t)Xj, ||6|, ^, pX| ]. 
The theory is then described by the Poisson, continuity and Euler equations, which are given 
by: 

= AnGp, (1) 
p + pvi,i = 0, (2) 

Vi + (^,^ = 0, (3) 

where G is the Newtonian gravitational constant. Here • refers to the convective derivative. 

We find from (^) that the Newtonian potential cp, up to a constant, takes the form (p = 
aij{t)XiXj. Thus equation (|l|) simplifies to 



an = AirGp. (4) 

We make the following decomposition 

Vij = -^OSij + (Tij + ujij, (5) 

with 

e = Vu, 

where the trace, 9 is the expansion, the trace-free symmetric part aij is the shear and the anti- 
symmetric part uJij is the rotation. Inserting these expressions for Vij, equations (|2|) and (^) 
become 



p + p9 = 0, (6) 



and 



e = 47rGp-^9^ + 2u}^ -2a^, (7) 
2 

LOi = --9uJi + aijLUj, (8) 
o 

2 2 1 

O'ij = —-^^ij — O'ik'ykj — i^ik^kj ~^ '^i'^ ~ ^ ~ '^ij + '^'^kk^ij , (9) 

where = ^aijaij and uj"^ = ^u>ijLOij = uiiUJi. 

Equations (Q) and along with a barotropic equation of state relating the density and 

pressure, form the Newtonian approximation to cosmology. The structure of these equations is 
as follows: There are eleven equations, one each arising from the Poisson, continuity and the 
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trace of the Euler equation, five from the symmetric trace-free part and a further three from the 
antisymmetric part. However, there are 16 unknowns: 6, p, tOi, aij and aij. Thus, the system is 
under-determined and it is not possible to solve uniquely for all the unknowns. 

Providing five arbitrary functions of time, for instance by specifying the shear (Tjj(t), for 
all time, we are able to solve for the five functions 0, p and w, the Newtonian potential (j) = 
aij{t)XiXj and the pressure p through the barotropic equation of state. Then, for a set of initial 
data (t(;j(xi, 0), p(io), i;^(io)) the set of homogeneous Newtonian cosmologies may be found since 
dH), (0) and can be uniquely solved for ^, ujij and p, and, an can be determined from @ and 
aij from (^). 

This is very different to the fully general relativistic case where there are nine possible 
symmetry groups. In general relativity the system of equations which describe cosmology has a 
well-posed Cauchy problem. In the Newtonian case, however, the system is under-determined 
- at least five functions of the set 9, uji and aij must be specified for all time. Thus, the 
Newtonian theory is not able to fully reproduce the full Bianchi types of general relativity in 
general. Arbitrary shear has no analogy in general relativity. However, once five (arbitrary) 
constraints have been imposed, the theory has an initial value formulation within the constraint 
of being homogeneous and isotropic and within the strange set up of allowing a given evolution 
of the shear. 

Let us explore how far the post-Newtonian approximation Q goes in overcoming this diffi- 
culty. The post-Newtonian theory is known to be well-posed in the general case of anisotropy 
and inhomogeneity |^. The extra potentials of the post-Newtonian theory act as integrability 
conditions for the Newtonian potential (j). We will see how these extra potentials help us in 
making the system determined in general. 

2.2 Post-Newtonian Approximations of Anisotropic Homogeneous Cosmolo- 
gies 



Following a scheme similar to that of Weinberg |12|, we adopt units in which the typical velocity 
has magnitude 1, i.e. /? ^, and assume a one parameter family of metrics g^j,y{x'^,c) for which 
there is a system of coordinates (x*^, x*) in which the components of the metric have the following 
asymptotic behaviour as c — > oo: 



3 I I /■|'^-^ 



9ij = Sij - 2(f)6ijC~'^ + aijC~'^ + a'ijC'^ + a-jC~^ (10) 

The usual Newtonian theory is obtained as the 0{c~'^) limit of ([lO|), while the Newtonian 
approximation is the 0(c~^) limit. Reformulating the field equations of the 0{c~^) limit as 
wavelike equations defines the post-Newtonian theory used here (see Q). 

After applying the harmonic gauge conditions 

<j> = -lCi,i, (11) 
Ci = 4>ij,j, (12) 

the field equations for the post-Newtonian theory are 

(p,kk = '^irGp + ^ {-4>jk,jk - J) , (13) 

Ci,kk = 16TTGpvi + ^ [ct)ijj - , (14) 
<i)ij - c^<Pij,kk = Kij + [16TTG{pviVj + 6ijp) - Jij] , (15) 
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where the a and aij of (10) are such that 



a 



26ija, 



with 



J = 



^Cj,j<P,i + 2C 



levrG 



ViP + pViV 



1 



"2 iCiCk,kj + CjCk,ki) 



- ^ijiQ<P,k(p,k + 327rG/)0), 

Cfc {Ci,jk ~l~ Cj,ik) ~l" '2CkCk,ij ~l~ Ck,iCk,j ~l~ Ci,kCj,k 
-2(j)^k i<Pki,j + '/'fci.i - 2(/>y,fc) - 16(/>(/),i0j + (t),i(t)kk,j 

+4>,j(t>kk,i — {4>ki,jk + (Afcj.ifc — </'ij,fcA: — (t>kk,ij) — 2(f)ki(t>,jk 

/I 1 1 

— 24>kj4>,ik + 2<j)^ij(j)kk — Sij [■:j^C,k,mCk,m + ^Cm.fcCfc.m + 2(CA:,fc) 

— CkCm,mk — '^4>,k4>km,m + ^4>,k4>mm,k — 12(/></>,fc0,fc 

— {2(l>km,mk — 2(f) mm, kk) 



+87rG 



2ptijT;j + 2p(20 + v^)viVj + 
1 , 3 , 1 



'+^(P,kP,k + ^, 



\kk 



This system forms a closed set which is consistent because the Bianchi identities are obtain- 
able from the field equations: 



P 1 + 



p [ 2vjVj + 2vjVkVk,j + ^Cj,j 



(16) 



and 



P{Vi + VijVj + + P,i 



1 



167rGc2 



-(j + Kjj)t;i + 



- Kijj - 20(Ji + (pjk,kij) - <P,i{J + 0jfc,jfc) 

Homogeneity can be provided by the following ansatze: 

(f) = aij{t)Xij + a{t), 

Ci = bijki{t)Xjki + bij{t)Xj, 

mn + Cijkl{t)Xki + Cij{t), 



(17) 



(18) 



where Xij = XiXj, Xijk = XiXjX^ etc. These are more general ansatze than in Q where 
isotropy has been ensured through similar constraints. As before we assume that the density 
and pressure are functions of time only, and are related through a barotropic equation of state. 
The velocity field decomposes as shown before in equation (|5|). 
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Substituting the ansatze (^) into the field equations for the post-Newtonian theory, Eqs. 
(13 to |l^), and comparing expressions of different powers in Xi, yields the set of equations for 
the anisotropic homogeneous post-Newtonian approximation: 



Qkk = 27rGp + c 



1 



Ckiki + SfrGpa 



(19) 



3 {2aikajk + Ckikiij) - 47rGp {Aaij + VkiVkj) = 0, 



(20) 



^bikkj = S-KGpVij + c 



Cikkj - 'ibkkaij - 2aikbkj + 2ajkbki + SvrG Vijp - -ph 



(21) 



2c, 



immjkl 



I kl I mm ™ J ™* 



kl 



+ 2a,mb^jkl - ^TtGp (^-b,,kl + ^y,^yrn[kfm{l] 



(22) 



'^Cijkk - IGaaij + 167rG(2pa + p)6ij + c ^ [-Cij + 2bikbjk + '^aijCkk - ^a-ikCjk - "^ajkCik 



-^ij [ ^bkkki + bkibki) + SttG (pcij + 6ij (-^ap + ^pckk 



0, 



(23) 



amkCLnk) + 32TTGp6ijar^ 

+8ttGp (VimVjn + VinVjm) 
2 



+C 



Cijmn ~\~ 16(2 ('iSijClkm^kn 2(lijyi(ljn 2(lin(i'jm) 



+ 2 [2b,nbkkrm + 2b,„bkk^m + 2b,^bkkjn + 2b,^bkk^n b,kb,kmn b.khkmn 
-\-2 {ciinCkkjm ~l~ ^jnCkkim ~l~ 0,imCkkjn ~l~ ^jm,C-kkin) 4 (dikCkjmn ~l~ (^jk^ihmn ^ij^kkmn) 
+ATTG(^2pCijmn + pVimVnj + pVinVrnj + pVjmVni + pVjnVmi 

^ij^mnkk 

= 0, (24) 



^ijmnrs 36^ | rnnr f ^kkj [7] | rnnr 



^ppi I s I ~^ "^^jfc I mn l ^j I rs | A; ~^ ^^jk | mn | rs | fc 



"64aj [^Oj |ir|°Prr| + ^^i [^'^fcfcj ["rarTI + ^'^j [mj'^kki \nr 
"^flj [^Cj-fcfc [H;^ + 4aijCfcfemnrs " (^ij ( 3 ^ rnnrs \ bkkll + 3b,.f^ [T^ hi j^jj - 4:8a^^^a^^akk 



-SttG 



|T]^j [7] + PCijmnr-s + 2pVi i^VjTT] 



( + 2/ja|^^VQr|fcV'fc|-j] 



0, 



(25) 
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where we have introduced the notation 



I mnr ^ ppj | s | ^^*"™''^PPi'' ~^ bifyirgbppjn + bimnsbppjr ~l~ binmrbppjs 

~^binr sbppjm ~\~ binmsbppjr ~l~ birmnbppj s ~l~ birmsbppjn 
~\~birnsbppjm ~l~ bismnbppjr ~l~ bigf^^bppjn ~l~ big^rb' 



'ppjm 



and similar, to indicate total symmetrisation in all indices that are surrounded by a box 
The harmonic gauge conditions (0^) give 







d 


= -^bkk, 






aij 




kj 




ijkl 





(26) 

and the Bianchi identities become 

L{t) + c-\M{t) + N,j{t)Xij) = (27) 

and 

OijX, + c-2(Qi,(t)X, + 5i,H(t)X,H) = 0, (28) 

where 

L{t) = p + pVkk, 

1 3 

M{t) = -4pa + -pbkk + pVkk - -^;^akkbu, 



Nij{t) = p (VikVjkVii + VikVjk + VjkVik + VikVkiVij + Vj-fc^fc/Fii + ^bkkij 

1 . 1 [7 7 

+ 2^ (^ifcVfcj + VjkVki — Sttij) + ^^^^ -aikbkiij + -ajkbkiii — '^^akkhiij 

Oij{t) = Vij + yjfcVfcj + 2aij, 

Qij(t) = ^l^ii (^4/9a + 4/9d - ^pbkk - paVkk - + 2Apaaij - pVij + + ^pkj 

P ] 

-{P + pa) VikVkj + -aij + 2paVkiVkj - pcikkj } 

~^ l&TrGp I ~ + Qbkkaij + Qbkkaij + 4aifc6fci + ^aikbkj - dGaaikajk 

—A&ackikiij + Sbijbkkii — 'ibikbiikj — Ghkbkuj + Gbkjbiku — Gbkibukj 
+^bkkbiiij — 9bkjbiiik — Qbkibkiij + Gbkibikij — 12aijCMkl + SaijCkkll — 8aikCkiij 

+^a.ikCiikj + 4:ajkCiiki + SukiCkUj — SakjCuki — 4afcfcq/jj| 



2 * J mm 



kl 



2Cr 



■immjkl 



V.^-^V mV 



[7] m ffc] m|T| 



.0 



+- 



48a i—id r— I + ISOmmb 



24a 



I k I mnn |T] 



+18a.|— 16 rrn + 18a.|— 16 rm + 4aim& rTTi + 4dimb 

" " " " " m\jkl\ n 



I j I mm I fc/ I i I j I mm | kl 

+48a, \-—\a m + 64a, \—\amma 



jki 



kl 



120a. I— ic 



0' 



kl 



[7] m|T] m|T| i [7]' 

48a;m.C j ^. I I I '^(^mmC-nnijkl ~\~ ^^(^imCnnmj kl ^^OjimCmmnj 



~\~ ^QO'mnCmnijkl 
kl 



+24a.|— ic 

« M mmnn 



kl 



-186 .|— 16 

mm J mn 



kl 



9b 
18h 



jk 



b [71 + 36j,fci6, 

m nnm \ I \ 



'mnni Omjkl 



'ijkl'^mmnn i ^^imnn^mjkl 
'2'^bmmnibnj kl 



+276 . |— 16 

mmi J nn 



kl 



+ 186 . |— 16 

mm J nm 



kl 



Note that setting L(t) = 0, i.e. setting higher order terms in to zero, one recovers the 
continuity equation of the Newtonian approximation. Similarly, 0{t) = yields the Euler 
equation of the Newtonian approximation. 

It can be shown that the time derivative of (1T9|) and the trace of equation (^T|) combined 
give 

L{t) + c-^M{t) = 0, (29) 
in accordance with the Bianchi identity ([l^). This essentially is the continuity equation with 



(30) 



c corrections. The time derivative of equation (^Oj) and the trace of equation ( |22D yield 

Nij{t)Xij = 



which is also consistent with the Bianchi identity (^). 

We can also recover the post-Newtonian Euler equation by combining the time derivative of 
( [2l| ) and the diagonal elements of (p^ to obtain 



The final part 



Oij{t)Xj + c-^Q,j{t)Xj = 0. 



SijkiXjki{t) = 0, 



(31) 



(32) 



can be obtained by combining the time derivative of (^) with equation (|25D. 

Thus we may completely define the anisotropic, inhomogeneous post-Newtonian cosmology 
with the set of equations; (||), (^, (|19|), (|2^), (g^), (|2|), (|5|), (p|, m / n) and In 
the case of the isotropic, homogeneous post-Newtonian cosmology, [0] it was shown that there 
existed relationships between the field equations rendering many of them redundant. This is 
special to the case of isotropy and due to the symmetries of the potentials. We know very 
little about the symmetries of the potentials in the more general situation of anisotropy. The 
harmonic gauge conditions provide a little assistance. For example: 

~ ~ ~ (11' (H 

191)7(1291) and 



Newtonian cosmology is contained within the equations (p9|), (plD, (|19|), (poD, ( |21 
(24, m 7^ n) and ([2^), as the special case when the terms go to zero. Equations ( 



is the only equation not to contribute to the Bianchi identities. 
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( |3l| ) are then the Newtonian theory with c ^ corrections. The corrections to the theory contain 
extra unknowns - the potentials Ci = hjkiXjki + hjXj and 4>ij = Cijkimn^kimn + CijkiXki + 
which are determined by the equations (^), (p|), (p^), (p^), (p4, ra^ ri) and (]2^ ). 

We win now attempt to use the equations (||), @), (111), d^), (H), (|5^, (H, m / n) 
and (H), to solve for the unknowns; aij(t), a{t), bijkiit), bij{t), Cijkimn{t), Cijki{t) and Cij{t). 

Let us start with equation (|2^) which has trace 



36, 



'Ikkl 



SnGpVkk + c ^ cikki - 'iakkhi + 8vrG 



pVkk - ^pbkk 



(33) 



With the aid of the harmonic gauge conditions (2C), it can be shown that this is just the time 
derivative of equation ([l9|). Thus this equation is redundant. The symmetric part of ([2l|) is 



(34) 



3 4 . 

2 {hkkj + bjkki) + -^akk^ij = SirGpaij 







{i 





2 ... 
- a 
3 



Apaij - p {bij + bji) - -paSij) 



where the harmonic gauge conditions and equation ( |19|) have been used. The antisymmetric 
part is given by 



S-nGpujij + - {bjkki — hkkj) 



+c 







{i 


bij bji 



'^aikbkj + 2ajkbki + 2-kG [Apuij - p {bij - bji)] > = 



(35) 



where once again the harmonic gauge conditions have been used for simplification where ever 
possible. These equations can be solved once we know aij to obtain information about ftjj^ and 
bijki- 

We now want to find aij. To do so we consider equations (pO|), (^) and (|25|). Equation (|2C 
has the trace 



akk = ^CLkim - 4vrG/9 ^^akk + ^O"^ + 2cr^ + 2^^! . 
providing us with an equation for akk- The traceless part is given by 



(36) 



Gaikttjk - CLij - A-kGp (^aij + ^O"^ + 2fT^ + 2u? + '^Ouij + aui^kj + (^kj^k^j = 0, (37) 



which yields aij. Summing over A; = Hn (22) gives an equation which is just the time derivative 
of (13). 



We still have information remaining in equation (22), however. The remainder of equation 
1221) is 



hjkl - 18«,g]'^| kl\m^ 



12a |— 16 . 

m J mi 



kl 



+ ^aimbmjki - SvrGp 



(38) 



Prom this equation the traceless part of bijki can be determined. 
Consider equation (|2^). The trace is given by 



2ckkU - IQaakk + 487rG(2/)a + p) 

+c"^ (-Cfcfc - bkibki + 4afcfcQ/ - SukiCki - 24d^ + AttG [5pckk - 3ap]) = 0. (39) 



^bkk of course being given by the harmonic gauge conditions. 
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It can be shown that, when combined with (jl^), this equation is just (^) again and thus is 
redundant. Now consider the symmetric piece 



16 I 1 

2cijfcfc - ^Cmmkk - IGaaij + —aauk^ij + c~^< - Cij + -Ckk^ij + '^hkhjk - ^hibkAj 



4 8 
+4aijCkk - -^akkCiidij - 4:aikCjk - '^ajkCik + -akiCkAj + SirGp 



= 0, 
(40) 

which can be used in equations (34) and (^) above to define Cij. There is no antisymmetric 



equation. 

We have seen that equations ( p^ ) and (|37| ) give aij{t). The harmonic gauge conditions can 
be used to determine bkk, and equations (|34|) and (^5[) give the trace- free part of bij We may 
then determine bijki in the fohowing manner: Use the harmonic gauge conditions to obtain 
bkkij (recah that there is symmetry in the last three indices). Then use Eq. (^) to obtain the 
remaining hijkis. The Qj's come from equation (|40|). To solve for Cijki use the harmonic gauge 
conditions to find Cikkj (where there is symmetry in the first two and second two indices) and 
Eq. (p^). Finally, Equation (p5|) then provides Cijkimn- 

Thus, equations (^), (|2T|), (p2|), ( |23|) a n d (|25| ) determine the unknowns aij, bijki, bij, Cijkimn-, 
Cijki and Cij, and the set of equations (|19|), (|29|), (|3l|), now contain only the eleven unknowns 9, 
iOij, aij, p and A(t). This means that the equations can be solved for all the unknowns uniquely. 

Therefore the post-Newtonian approximation provides a well-posed, closed, complete system. 
Thus, the post-Newtonian approximation produces a set of cosmological equations which are 
more in the spirit of the Bianchi types of general relativity. 

3 Newtonian Theory and the Post-Newtonian Approximation 
of Rotation-free Anisotropic Homogeneous Cosmologies 

3.1 Rotation- free Anisotropic Homogeneous Newtonian Cosmology 

We first like to examine the Newtonian approximation of cosmology with vanishing rotation and 
diagonal shear, aij = diag(c7ii, o"22, — cii —(^22)- The Newtonian approximation (Eqs. (Q) and 
(|6|) to (|9|)) is then given by 

an = iirGp, (41) 
p + pB = 0, (42) 

e = ATrGp-\e^ -2a\ (43) 

2 2 1 

<jij = aij - -9aij - aikakj + 6ij{-a - -akk)- (44) 



For i ^ j, equation (|3) yields 



0, i / j. (45) 



The diagonal elements of the shear obey the differential equation (44, i = j, no summation 
over i) 

o'ii = an — -Ban — 'Ek'^ik'^ki + — -^kakk- (46) 
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Summing over i just yields the definition of o"^ 

2cr^ = (Tijaji = 2 (all + ^22 + ^^110-22) • (47) 

Thus, ( |4^ ) provides at most two independent equations. Introducing 9 = 3-^ the set of equations 
can be written as 

R A 2 

- = -TrGpoR"^ - -{ali + al2 + an(T22), (48) 

P = PoR-^ (49) 

R \ 2 4 

On = 0-11 + 2— cTii + -all - 0(^22 + 0-11^^22) + -ttGpqR'^, (50) 
rt o o o 

7? 1 2 4 

022 = 0-22 + 2—0-22 + -o"i2 - o('^ii + '^11^22) + -irGpoR'^, (51) 

033 = 47rGpo^"^ - an - 022, (52) 
where po is a constant. 

Consider the structure of these five equations: If we provide the functions an and a22 for 
all time, then the theory has a well-posed initial value problem for the variables R{t), p{t), 
aii{t), 022 (t) and 033 (t). Equation ( ^8| ) is the Raychaudhuri equation (obtained from (|43|)), 
and its solution provides R{t). From Eq. (|49| ) (which comes from (^2|)) the density /9(t) can 
be extracted which, once an equation of state is provided, gives p{t). Equations (|50|) and (|5] 



(from (l44|)) give an and 022, and from (|52D (or (|4l|) ), which is the Poisson equation, 033 can be 
extracted. 

Thus, for a homogeneous anisotropic cosmology in Newtonian theory, there is an initial value 
problem provided we supply two functions for all time. We have seen that the post-Newtonian 
theory is able to give as many equations as unknowns. So there is no need to specify any functions 
of time. Actually, there is no need to (arbitrarily) specify any of the unknown functions in order 
to obtain an initial value problem. We will revisit the Newtonian theory shortly, when we will 
consider the Bianchi I metric as a specific example of a homogeneous anisotropic cosmology with 
vanishing rotation. 



3.2 The Homogeneous Bianchi I Universe 
3.2.1 The Metric 

Now we would like to study Newtonian theory and the post-Newtonian approximation as applied 
to a specific example of an anisotropic homogeneous cosmology. To this end we consider the 
homogeneous Bianchi I models - i.e. the Heckmann-Schiicking solutions of general relativity.^ 
The metric of the general Bianchi I universe has the form 



ds'^ = —dxQ + RikRjkdxidxj, where xq = ct with Rij = diag(i?ii, i?22> -^33)- 

To proceed further, we need to write the Bianchi I metric into a form off which the potentials 
(f), Q and (pij can be read. To do so, we consider the following transformation 

Xo = Tc + TC~^ + t' , 
*Not to be confused with the Heckmann-Schiicking cosmology of Newtonian theory. 
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' -4 



(53) 



with 



r = A{t) + A^j{t)Xij, 
t' = Bij{t)Xij + Bijtti{t)Xijki, 
Xi = Cij{t)Xj + Cijki{t)Xjki, 



Xi 



Dijki{t)Xjki + D, 



ijklmn 



{t)X 



jklmni 



where Xij = XiXj and similar for Xij^ etc. Throughout the remainder of the paper we assume: 
j4 ^ to ensure convergence of the expansion in c~^. A{t), Aij{t), Bij{t), Bijki{t), Cij{t), 
Cijki{t), Dijki{t) and Dij^imnii) arbitrary functions of time. In these coordinates the metric 
reads 



l + c" 



+c 

-2i ( R^R-^ 

+cdTdXi "-^ 

+c 



A + AijXij 



2BijXij — IBiikiX, 



ijkl^ijkl 



Xi 



2 Rim, Cm. iki. X', 



■i'm'^mjkl^ijkl 



2RjiCiiXij 



Xi 



2Aki [R'R-')^^Xijki\+0{c-^) 



-4:A,jXj - 2 ( RR- 



SAAijXj — SAijAkiXjki — ABijXj — BijkiXjki + 4Aj; (^R^R ^^^^ Xjki 



-\-2RjiCilXj + 2RiiClj — GRjmCmikl^jkl + '^RimCmjkl^jkl 



+2A[RR),, Xj + 2Au [RR] , X^ki] + 0{&) 



+dXidXj 



+0 



6ij + c 



^AikAjiXki + 4:Ajk [RR] ,Xki + 2RikCjk + QRimCmjkiXki 



SAAikAjiXkl - SAmnAikAjiXklmn — ^BjkAuXkl - IQBiklmAjnXklr 



~\~^RimBfnjklXkl + ^ORipDpjklm.n-^klmn ~l~ RranC imC jn ~l~ ^RipRjqCpiklCqjmn-^klmn 



+'^RmnCimCnjklXkl + "^RmkRknC jnCraiklXkl + 4 ( i? j AimAjnXkl 



~^^AjiRkmCmi^kl ~l~ ^'^Aj^RkpCpHmXklmn ~1" '^Ajm,RilCklXkm ~\~ '^Aj^RpiCpklm-^khi 
+4i {RR'^)^^ AjiXkl + 4 {RR~^)^^ AranAjlXklmn + 4 (^^"')^^ BjiXkl 



+8 [RR ^ ) BiimnX, 



ik 



'jlmn-^klmn 



0{c 



(54) 



It can be seen from (|lO| ) that there are no terms of order c ^, thus 

1 



RR 



-1 



Also, from (|lO|) we see that the c ^ terms in goo and gij must be identical. Hence, 



Cijki 



RR 



-2 



RR 



-1 



jk 



+ ( RR-^ 



ik 



RR 



-1 



+ \r^I {rr-^ 



kl 
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We may now read off the potentials </>, Q and 4>ij from the metric (^). With the help of the 
harmonic gauge conditions (0) and (|l^ ) we are able to write down the functions a(i), aij{t) 
etc. which appear in the ansatze (|l8|): 

a{t) = A, 

1 



bijki{t) = -SBijki + — ( 3(5 



36 



kl 



iiR-^ 



-5 



-\-3RifjiD.i 
+2Bji (^RR-^ 



J ji 2 \ ni 2 V J jk 2 \ Jik 

^mjkl ~l~ 3RjrnDmikl ~\~ 2Bik 

RR-') -^A(rr-') (rr-') -^A(rr-') (rr-') 

J ik 4 V Jik\ J jl 4 V J il \ J jk 



il 



V /i/c 4 V /i/fc V /ji 4 



4 

4 V 



/ jk 



(rr-^ 



, --A( RR- 
'ii 4 V 

1 



1. 



mn I 



kl 



I mn I 



+25, 



(rr- 



kl 



, , - - (rr-^^ 

mn 3 V 



kl 



(rrr~ 



I mn I 



— (5jj ^ {^R R ^ I I + ^RioDojklmn + ^RjoDoiklmn 

+8B 



1 
6 

25 



+ -<5 r-^ (RR-^)r-i (rRR-^), , --S r-^ ( RR^^)r-i (r R"^), , 

6 i\k\\ /[jljv / [mn] 6 i\k\\ J [mn\ 

+ (^^")^ + (^^"O.a (^^"^ 

4 (^^").m - 1 (^^").m (^^ 



.0 



I mn I 



R^R-^) 



mn I 



I mn I 



- 2 



(55) 
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The Heckmann-Schiicking solution has Rij = diag (i?ii, R22, R33) such that 

Ru = {t-to)P'{t-hr, (56) 

R22 = (^-^o^(^-^lr^ (57) 

R33 = {t-tont-tiy^ (58) 

with 

EiPi = ^ipj = 1, (59) 

where 

2 

qi = ^-Pi (60) 

such that 

Sj^i = Sj^f = 1. (61) 
These conditions can be fulfilled by the following parameterisation for pi to P3 and qi to q^: 



Pi 


:= u, 




l-u 1 


P2 






l-u 1 


P3 


2 2 




2 


qi 


:= u, 

3 




1 + 3n 




6 




1 + 3n 


q3 


^= 6 + 



Vl + 2n-3n2, 



-Vl + 2u-3n2, 



(62) 



3.2.2 The Heckmann-Schiicking Solution of General Relativity and Newtonian 
Theory 



The Bianchi I model is rotation- free and so the Newtonian theory is just given by equations (45) 
to (^). From (|^) we see that = — Define R{t) such that 

g _ _ {R11R22R33)' _ 2t — tp — ti ^g^^ 
R R11R22R33 {t — to){t — ti) 

The parameters to be determined are n, to, ti and po, and the variables are p, an, a22- 
First consider equations (|50|) and ( [5l| ) which using the parametrisation ( |62[) become: 



2 2t-to-ti 1 2 2, 2 X 4 ^ 1 



3 {t-tQ){t-ti) 3 3' ' 3 '^(t-to)(i-ii) 

u(n-l) _^ (2-3n)(l + 3ti) ^(2-3^) 



2{t-toY lS{t-tif 3(t - to){t - h) 
and 

2 2t-to-ti 1 2 2, 2 N 4 ^ 1 

0"22 + -<7227- TTTL TT + 0*^22 ~ ol'^ll + '''ll'^22j + o^^/''; ~ 



3 {t-to){t-ti) 3 3' ' 3 {t-to){t-ti) 
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(l - n + Vl + 2n - 3ti2) (l + n - Vl + 2n - 3ti2) 
8(t-to)2 

(l + 3it - 3Vl + 2u - 3^2) (5 - M + 3Vl + 2n - 3ti2 
^ 72(t-ti)2 

(1 - n + Vl + 2n - 3^2) (1 + 3-u - 3Vl + 2u- 3n2) 
12(f-to)(t-ti) ■ 

Once the parameters u, to and ti are known we are able to solve these two equations for an 
and (T22-Q To proceed further, consider Eq. (48) which is a Raychaudhuri-like equation and may 
be written as 

4 1 2,22 ^ 1 2 2t-to-ti 

o'^POTl — 7171 — TT ~ o (^11 + '^22 + f'"ll'^22 



{t-to){t-ti) 3' 3((t_to)(t_t^))f 9^(^t-to){t-ti))l' 

If cJii and (T22 are known then this equation relates the two constants to and ti. The Poisson 
equation (^2|) simply defines po, 

12^- 

Finally, we can write down ( |49| ) as 

1 



(i-to)(t-ii)^ 



which determines p{t). Therefore the shear variables, an and (T22, and the density p are now 
expressed in terms of two parameters with the initial condition p{t = 0) specifying one of them. 

Thus, the Newtonian theory in the homogeneous, anisotropic rotation-free case can be solved, 
provided we arbitrarily fix the parameter u for all time. However, varying the equation of state 
will have no effect on the solutions for the density, since the pressure in no way enters the 
dynamics, as was the case in Q with the FRW cosmology. 

3.2.3 The Post-Newtonian Approximation of the Homogeneous Bianchi I Universe 

The post-Newtonian theory for the homogeneous anisotropic rotation-free case is given by equa- 
tions (1^), (H), dH), dH), (H, m / n) and (|3|) which may be solved for the unknowns ciij, 



bij, bijki, Cij, Cijki and Cijkimn as we have shown in section 2. The set of equations (p^, (2£) 
and ( |3l| ) then determine the eight unknowns 6, aij, p and A{t) (or a{t)). Vanishing rotation 



imposes extra symmetry on the equations. Therefore many become redundant and there are a 
great many simplifications in what does remain. Also we can use the potentials (|55| ) for further 
simplification. Clearly aij = for i 7^ j, since aij = — with Rij = diag{Rii, R22, -Rss)- 
Using the parametrisation (^) the Poisson-like equation (|T^ ) becomes 



^ evrGp + 2^T^Gp + ^ ^ I , (64) 



{t-to){t-ti) 



which can be used to determine the potential A{t). 

The Friedmann-like equation (p9|), as usual, provides p{t): 



2t-to-ti , A . .. 2t-to-ti , A 1 1 . ...^ 

^ + ^t-^o)(t-t,)+^ - - 2pA + p^^_ ^^^^^ _ + ^ _ _^A = 0. (65) 



The condition an + an + 0-33 = determines (T33. 
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The Bianchi identity (^Tj) can be decomposed into a trace part, symmetric trace-less part 
and an antisymmetric part. The trace yields the Ray chaudhuri- like equation 

6{t - to){t - h) - 2{2t - to - tif 



{t-tont-hy 



+ 



-c-2< 



A 
9" 



14 



I24:{2t - to - tif 

(i-to)2(t-ii)2 " {t-to){t-ti) 



22 ■■ 2t — to — ti 



3 it-to){t-ti) 



123 ... 2 



/63 l {2t-to-tif ^ 
™"ll6 2{t-to){t-ti) 



5.3 

"o^ + o^kkllit - to){t - ti) 
6 o 



9, 



+ CkkllC^t — to — ti) — -bkl {bklmm + blkmm) 



-96 A 



+ 



u{3u - 2) 



u{u-l) ^ (2 -3u)(3n + 1) 



l8{t-toy 



u 



3{t-to){t-ti) 2{t-tiY 

1 - ^/l + 2u- ?,uA (-U - 1 + Vl + 2n - 2,v? 



4{t-to){t-ti) 



+ 



l-u + Vl + 2u-3u'^) (u + 1- Vl + 2u-3u'^ 



u 



8{t-tif 

(Sm + 1 - 3Vl + 2u- 3n2) (su - 5 - 3^1 + 2u - Su"^ 
72{t - toy 

1 + Vl + 2ti - 3n2) (3u + l + 3Vl + 2u - 3^2) 
I2{t - to){t - ti) 

(l-u- Vl + 2u- 3n2) (u + 1- Vl + 2u- 3n2 



+ 



8(t-ti)2 

(Sn + 1 + 3Vl + 2u- 3n2) (^3u - 5 + 3Vl + 2n - 3^2 



+ 24^- 



72(t - toy 

Qu (3n - 2) + Mu-m-to) ^ {3u - 2)(3n + l){t - ti) 



+ 



t — ti t — to 

6 (1 - + Vl + 2u - 3^2) (1 + 3u - 3^1 + 2^ - 3it2) 

9 (1 - u + Vl + 2u - 3^2) (1 + u - Vl + 2u - 3n2) {t - to) 

t - ti 

1 + 3n - 3Vl + 2u- 3ii2) (-5 + 3n - 3^1 + 2n - 3^2) (i - ti)' 



H ^ss 

96 



+ 



6 (1 - u - Vl + 2u- 3^2) (1 + 3ti + 3v/l + 2u - 3?/2) 

9 (1 - u - Vl + 2u - 3n2) (1 + u + Vl + 2u - 3n2) (t - to) 

t - ti 

(1 + 3n + 3V1 + 2u- 3ii2) (-5 + 3« + 3\/l + 2u - 3u^^ {t - ti) 

t^o 



21 
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1^ 

D 



3u(2-3n)(2t-to-ii) 9u{u - l){t - to) (2 - 3u){3u + l){t - h) 



+ 



H 622 

48 



12 



6 (-1 + n - Vl + 2n - 3u^) (l + 3u - 3Vl + 2u - 3^2) (2t - to - ^i) 

(t-to)(i-ii) 

(-1 - u + Vl + 2?i - 3?x2) (1 + u - Vl + 2u - 3n2) (t - to) 
(1 + 3m - 3Vl + 2ti - 3^2) (-5 + ti - 3Vl + 2u - 3^2) {t - ti) 



{t-tof 



1 



+ 777^33 



6 (-1 + n + Vl + 2m - 3n2) (l + 3m + 3Vl + 2m - 3m2) (2t - to - ^i) 



+ 



{t-to){t-ti) 

12 (1 - M - Vl + 2m - 3m2) (1 + m + Vl + 2m - 3m2) [t - to) 
(1 + 3m + 3Vl + 2m - 3m2) (-5 + m + 3Vl + 2m - 3m2) (t - ti) 



which may be used to determine u. 

The traceless symmetric part is given by 



0, 



(66) 



+ CTjj 



2 2t — to — t\ 



3 {t-to){t-ti) 



+ <7,-i — -c 6ij + 2a 



9it - to){t - ti) 



Si 



T2'''' 



IO6A - h^A- ^— + 3(t - to)(t - t\)ckm + 3 (2t - to - ti) Cfc^/ 



(t-to)(t-ti) 
A + 9A(2t-to-ti) 



70-1 4 • 

+— ^ajj + -(t - to)(t - ti)ckkii {23aij + 2&ij) + -A&ij 
00 o 



1 



27(t-to)(t-ti) 
+it-to){t-ti) 



Sii - - { bij + bji ) - —CkkiiS 



9 



23 



72 



. ^ . . 7. 

3 ^ Ojj + 9Adij + - {ciikbkj + ajkbki) — -akibkAj 



-SGauikajk + l2aakiakiSij + - {ajkhk + cukhjk) 



(67) 



— 7 {bikbkjii + bjkbkiii + "ihkibkuj + "ibkibikij — bkibkimmSij bklblkmmSij) 
3 

+ 2 {(^ikClljk + O-jkCllik) — akmCllkmSij + ^aklCklij 

and may be used to find cth, 022 and (T33 = — fxn — 022- 
The antisymmetric piece of Eq. (^l|) is 

^{^ikkj — ajkbki) + 2{djkbik — dikbjk — ajkhk + aikbjk) 

= '^ibikbkjii — bjkbkiii — bkjbikii + bkibjku). (68) 

These are three equations. However, with the aid of the harmonic gauge conditions it can be 
shown that only two of the equations are independent. These can be used to determine to and 
ti. 
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Thus, the Heckmann-Schiicking solutions of general relativity for the post-Newtonian ap- 
proximation are defined by the eight variables /?, i?22) R33, cth, 0-22 o"33, which are 
completely determined by the equations ( |6^ to (|6^). 

A second improvement on the Newtonian theory is that here pressure enters into the dynam- 
ics through the Raychaudhuri equation. This means that varying the equation of state varies 
the solutions for the density p, the cosmic scale factor R{t) and the shear cr{t). As in the FRW 
cosmology the function of time A{t) is what incorporates the pressure into the system. 

4 Newtonian Theory and the Post-Newtonian Approximation 
of Shear-Free Anisotropic Homogeneous Cosmologies 

4.1 The Newtonian Shear-Free Anisotropic Homogeneous Universe 

The Heckmann-Schiicking solution of section 2.1 is the Newtonian approximation of an 
anisotropic homogeneous cosmology. The cosmology is given by equations (Q), (^) and (^) 
to (P) with vanishing shear. The equations reduce to 



an = AnGp, 

P = PoR-^ (70) 
iOi = ujioR~^, (71) 
R _ AttGpo 2i4 

R~ 3R^ ^ R^- ^ ' 

Recall that with the five shear functions arbitrarily set to zero we have a determined system. 
The final equation (^) may be integrated out to give the Heckmann-Schiicking solution 

i?^ = ^-^-4, (73) 

where e = it 1 or 0, Tq an arbitrary constant. When Wjo = we have a shear- free and rotation- 
free cosmology which becomes isotropic. This is just FRW, which is the most general possible 



solution for homogeneity and isotropy. In this case equation (73) can be identified with the 
Raychaudhuri equation of the FRW cosmology. 

A theorem of Ellis [|l| states that in the case of shear-free dust either the expansion or the 
rotation must vanish. Setting ^ = 0, we obtain the following solution for R{t) 

SirGpo ^ ^ 

Hence i? is a constant, and as long as ujq is non-zero, there is no possibility of there being a 
singularity for R{t). The solution will have minimum and maximum values for R{t) but may 
never be zero, pO[ . 

Next we explore the analogous case in the post-Newtonian approximation and see how far it 
comes in overcoming this difficulty. 

4.2 The Post-Newtonian Approximation of an Anisotropic Homogeneous 
Shear-free Universe 

The post-Newtonian approximation for the shear-free case is given by the field equations (|20|), 
([2l|), (p^), (p^), ( p4| , m ^ n) and (p3|), along with a Poisson-like equation (|l^) and the Bianchi 

®And, p{t) which is defined through a barotropic equation of state. 
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identities ([29|), (pi]), with = 0. The equations (|T9|), ( |29| ) and (^Tj) are eleven equations, and 
the variables are reduced to p (which may be determined from (p9|)), (which may be determined 
from the trace of (|3ll)), oJij (which may be determined from the antisymmetric part of ([3l|)) and 
a{t) (which may be determined from the Poisson-like equation (^)). Unlike in the rotation-free 
case where the aij can be determined from the additional field equations here we will use what 
remains of the symmetric part of (^) to determine the for i ^ j. The unknowns hijkh 
Cij, Cijki and Cijkimn that appear in the corrections t o ([T9| ), ( |29| ) and (|3ll ) are determined 
from the additional field equations (20), (^, (22), (23), (|24 m / n) and (pH)- 

The reduction in unknowns must be accompanied by redundancies otherwise the system 



will be overdetermined. We saw earlier that equations (36) and ( p7| ) yield Uij. Somewhat 
tedious calculations show that the time derivative of the symmetric part of (^) and (^) give 
the symmetric part of ( |3l| ) in this special case of vanishing shear. Thus, we may ignore the 
symmetric part of (|3l]). 

The Poisson-like equation (^) is 

Ukk = 2-kGp + c~'^{-^Ckiki + SirGpa), (75) 

which defines a{t). Defining 6 = 3^j^ we may write the continuity equation ( |29| ) which allows 
us to solve for p as^ 



= ^P + 7c"^ (76) 



1 

(i?lli?22-R33)^ 



where 7 is a solution of the differential equation 



i\ ..\ / _ _ _ 1 

3 



7 + 47^^ + A 



^ 1 
(^11-^22-^33)^ (-^11-^22-^33)^ 

_(i?lli?22fi33)^)" ((i?lli?22i?33)^)' ((i?lli?22i?33)^^ 

+8d 2 i — — 0- (77) 



(i?lli?22i?33)3 {RnR22R33y^ 

The antisymmetric part of (|3l|) becomes 



+ '^(^^ij + ""^ ^ (^^i (^~2pa - 2pa + 2pa^ + ^pbkk^ - ^P (cikjk - Cjkik) 



4 ^ 1 
+2pau>ij — —paOiOij + 



3^ levrG 



-OuJij - -ujij (d +Adakk + 4aafcfc) 



4/2 1 1 \ 

-- ( + -6'a;jj + -LOikUJkj - -^^jk^kij (a + ^aa^fc) - 2Auijakkd - 4:dkjbki 

+4dfci6fcj — Aakjhki + ^akibkj — l^bikhikj + l^bjkbuki — 12aikCiikj 
+12ajkCiiki — 'iuikCkiji + ^ajkCkiii — -iajkCuki + '^aikCjiki + SaikCukj 
—8aikCijki — SaikCikij + SaikCjku]) = 0, (78) 



which may be solved for uJij. Finally, the trace of ( |3l| ) gives a Raychaudhuri-like equation which 



IS 



R ^ / 9\ 9 9 / R\ R 95 

3- = 47rG(p + 3pc~2j+2a;2 + c-2 -ii4q _ + 20d- + Saa^fc - — a 



'^Which is the same as in the rotation-free case. 
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--aaklttkl + akl {2Ckmlm — 8Cmmkl — GCmifcm) + Hi 



(79) 



Thus, the case of an anisotropic homogeneous shear- free cosmology is essentially defined by 
a Poisson-like equation (^) , a Friedmann-like equation ( |76D , a Raychaudhuri-like equation (|7^) , 
and has rotation given by (|7^). Its variables are p{t), R{t), u>ij{t), aij{t) and a{t). 

In contrast to Newtonian theory, in post-Newtonian theory the Raychaudhuri equation has 
the pressure entering the dynamics and so gives rise to a variety of possible solutions for p{t) and 
R{t). Before the only solution was that of dust. It is through the function a{t), that the pressure 
enters into the dynamics. Without such a term the c~^-corrections would vanish and we would 
be left with a system where pressure is not dynamic. There is however, a more serious reason 
for keeping these higher order terms and it relates to the Ellis theorem: The Raychaudhuri-like 
equation ([79|), allows for solutions with R{t) = 0. Due to the pressure becoming dynamic, with 
a variation in equation of state the solutions for the density vary. Thus equations (^) and (71) 



no longer hold. Solutions different to dust are possible, and the rotation is no longer just a 
constant. Hence the Raychaudhuri equation (|79D will no longer yield the result (|73). 



5 Conclusion 

In this paper we have derived the post-Newtonian approximation for anisotropic homogeneous 
cosmologies. In contrast to the Newtonian approximation the equations are well-posed. The 
cosmological equations of the post-Newtonian approximation are much more in the spirit of the 
Bianchi types of general relativity. 

We considered a particular example of the Bianchi identities - the Bianchi I cosmology. In 
the Newtonian theory assuming anisotropy and homogeneity leads to the well known Heckmann- 
Schiicking solutions. Howsoever, even with homogeneity, there are still not enough equations 
to solve for all the unknowns and some need to be supplied for all time. The post-Newtonian 
theory is able to overcome this problem and also allows the pressure to enter into the dynamics 
of the theory. Therefore the full set of possible solutions for R{t) as outlined in [Q], can be 
reproduced. 

The result of pressure entering into the dynamics allows for more than just a matter domi- 
nated universe. In this way more general solutions are obtained than in the Newtonian case. In 
the shear-free case with rotation the Raychaudhuri equation gives rise to a singularity. In this 
way the Ellis theorem does not lead to contradictions, and there is hope that the solutions of 
the post-Newtonian approximation may all have general relativistic counterparts. In the very 
least we no longer have solutions to which there can be no general relativistic analogues. 

It seems reasonable then to replace the Newtonian theory with the post-Newtonian approx- 
imation when considering homogeneous cosmologies, whether they be isotropic or not. 
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